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Abstract 
^ ■ 

\ Exact non-static spherically symmetric solutions of the Einstein equations 

qh| for a null fluid source with pressure P and density p related by P = kp"" 

^j^l are given. The a = 1 metrics are asymptotically flat for 1/2 < A; < 1 and 

cosmological for < < 1/2. The k = 1 metric is the known charged Vaidya 

X 

?H ' solution. In general the metrics have multiple apparent horizons. In the long 

time limit, the asymptotically flat metrics are hairy black hole solutions that 
'fall between' the Schwarzschild and Reissner-Nordstrom metrics. 
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The problem of gravitational collapse in general relativity is of much interest. One 
would like to know whether, and under what initial conditions, gravitational collapse re- 
sults in black hole formation. In particular, one would like to know if there are physical 
collapse solutions that lead to naked singularities. If found, such solutions would be counter 
examples of the cosmic censorship hypothesis, which states that curvature singularities in 
asymptotically fiat spacetimes are always shrouded by event horizons. 

Since the general problem appears intractable due to the complexity of the full Ein- 
stein equations, metrics with special symmetries are used to construct gravitational collapse 
models. One such case is the two-dimensional reduction of general relativity obtained by 
imposing spherical symmetry. Even with this reduction however, there are very few inho- 
mogeneous non-static exact solutions known. One example is the Vaidya metric It 
describes the collapse of pressureless null dust and is asympotically fiat. For the minimally 
coupled massless scalar field, the only inhomogeneous non-static examples are cosmological 
solutions which have naked singularities [0,0. For perfect fluids with equation of state re- 
lating pressure P and energy density p given by P = kp, self-similar solutions (which are 
not asymptotically fiat) have been studied 0. 

Recently there have been a number of numerical studies of the spherically symmetric 
collapse problem [|6|-[T0|1. The initial data in these works is a compact ingoing pulse. It was 



found that a black hole forms if any of the parameters c in the initial data are above certain 
critical values c*, otherwise the pulse scatters back out to infinity leaving fiat space. In 
particular, the black hole mass formula Mbh = k{c — c^,)'^^ was discovered in Ref. 0. This 
shows that black holes initially form with zero mass - that is, no mass gap is found. 

There have been a number of attempts aimed at an analytical understanding these results 



Tl|-fi8[. One approach is to look for exact collapse solutions and see what black hole mass 
formulas may be extracted from them. 

In this paper we give exact inhomogeneous and non-static spherically symmetric solutions 
of the Einstein equations for a collapsing null fluid. The equation of state of the fluid is 
P = kp"-. As we will see, while these collapse solutions do not have direct relevance for the 



critical behaviour mentioned above, they do have a number of interesting features, including 
hair on black holes. 

An inverted approach is used to find the solutions. First the stress-energy tensor is 
determined from the metric. Then the equation of state and the dominant energy condition 
are imposed on its eigenvalues. This leads to an equation for the metric function, which is 
easily solved. The precise form of the stress-energy tensor is then displayed. Some interesting 
properties of the solutions are then discussed. 

For the general spherically symmetric metric 

ds^ = -e2V'(-.^)F(r, v) dv^ + 2e'^("'") dvdr + r^dfl'^, (1) 

where < r < oo is the proper radial coordinate, —oo < f < oo is an advanced time 
coordinate, and dfl'^ is the metric on the unit two sphere, the Einstein equations Gab = ^nTab 
give 

m'^^ = -AnrX', (2) 
rn^^= 4v^rX^ (3) 
ip' = AurTrr. (4) 

The mass function m{r,v) is defined by F{r,v) = 1 — 2m{r,v)/r, and is measure of the 
mass contained within radius r. We will consider the special case ip{r,v) = 0, which means 
from §) that T„ = 0. 

The stress-energy tensor derived from the above metric may be diagonalized to give the 
energy density and the principal pressures. The eigenvalue problem is Tj'Uh = XUa- The 
9-(f) part of T^*, which is determined from Ggg and G^j,^ = sin^ 9 Ggg, is already diagonal 
with pressure eigenvalues 

F . r/ = T/ = (6) 

Since T^r = (by choice), we have T^^ = T/ and T^'" = 0. Therefore the v-r part of the 
matrix to be diagonalized is 



V T/y 



(6) 



This has one eigenvalue A which gives the energy density p, namely 



m' 



-P-A = T/ = -— . (7) 

The corresponding eigenvector is Va = (1,0,0,0) (in the coordinates {v,r, 9, (/))), and is 



lightlike. Therefore the stress-energy tensor, (which follows from ip = 0), is oi Type 11 ||T9 
Its non-vanishing components are 

2m m 

T^^, = p {1 ) + - — -, T^r = -p, (8) 

Tee = P gee-, ^</-0 = P 9H-: 

These components may be succinctly written using the two linearly independent future 
pointing lightlike vectors Va = (1, 0, 0, 0) and Wa = {,F/2, —1, 0, 0) as 



VaVb + p{r, V) {VaWb + VbWa) + P{r, v) {gab + VaWb + VbWa) • (9) 



(This tensor should be compared with the perfect fluid one Tab = P UaUb + P { Qab + UaUb ), 
where Ua is timelike.) 

The stress-energy tensor (|^) has support along both the two future pointing lightlike 
vectors Va and Wa, and as we will see below, it is of precisely the form which gives the 
charged Vaidya solution 0. For P = p = 0, (^ reduces to the stress-energy tensor which 
gives the uncharged Vaidya metric. We also note that while TabW°'w^ = rh/Airr'^, TabV°'v^ = 0. 
Therefore there is energy flux only along one of the null directions. 

A static observer with 4- velocity S*" = (1/-\/F, 0, 0, 0) and a rotating observer with 



4-velocity i?" = (^2/^, 0, 0, 1/r sin 6^) see respectively the energy densities TabS^S^ = 
Th/AirFr^ + p and, TabR'R^ = 2{m/ATrFr^ + p) + P 

The stress-energy tensor (|^) satisfies the dominant energy condition if the following three 
conditions are met : 
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P>0, p>P, and TabW^w^ > 0, (10) 

The first two of tliese imply tliat m' > and m" < 0. Tlie former just says that the mass 
function either increases with r or is constant, which is a natural physical requirement on 
it. 

To satisfy the first two of the dominant energy conditions (0), we impose the equation 
of state P = kp"-, with k < 1 and a < 1 . The a < 1 case will be discussed later. For a = 1 
this gives the equation 



m" , m' 



(11) 



for the mass function, which is easily integrated to give 

fiv)-giv)/[{2k-l)r'''''], k^\ 
m[r, V) = > [12) 



f{v) + g{v) Inr , k 



where f{v) and g{v) are arbitrary functions (which are restricted only by the energy condi- 
tions). Therefore we have explicitly that 

^' = *i^ = *P- (13) 

Therefore we must have g{v) > for positive pressure and energy density. 

The last requirement in (|T0[) for the dominant energy condition leads, for /c 7^ 1/2 to 

"-^W- (2.-l)U-' ""- 
Physically this means that the matter within a radius r increases with time, which corre- 
sponds to an implosion. This condition is most easily satisfied if / > 0, and either g > 
and k < 1/2, 01 g < and k > 1/2. 

In summary, for the dominant energy condition, we must have g{v) > and, either g > 
for k < 1/2 OT g < foT k > 1/2. For the weak or strong energy conditions, (which are 
equivalent for Type II stress-energy tensors), we only need p > 0, P > and Tahw'^w^ > 0, 
but not p > P. Therefore for the latter energy conditions we can have k > 1 as well. 



For k = 1/2, neither the weak nor dominant energy conditions can be satisfied for all 
r because rh = f{v) + g{v)lnr, which always becomes negative for sufficiently small r. 
Therefore we will not consider this case further. 

In summary, we have shown that the metric 

^ _( 1 _ '^M. + )dv^ + 2 dvdr + rHn-" (15) 

r [2k — 1) r""* 

is a solution of the Einstein equations for the null fluid stress-energy tensor (H) with P = kp, 
where P is given by ([131) . 

There are two special cases of this solution which are already known. One is the Vaidya 
metric [|l|, which arises for g{v) = (vanishing p and P). Then the only non- vanishing 
component of the stress-energy tensor (^) is T^^ = m/Aiir'^ = /(f)/47rr^. The other is the 
charged Vaidya metric, where the charge depends on 0. This arises when A; = 1 in ([T5|). 

We therefore see that the stress-energy tensor we have determined by imposing the 
equation of state P = kp is a one parameter [k > 0) generalization of the stress-energy 
tensor which gives the charged Vaidya metric. The corresponding metric must therefore 
depend on k as in equation (P^). (We note that the parameter k in the metric or the 
stress-energy tensor cannot be eliminated by a coordinate transformation because it is the 
constant of proportionality between the eigenvalues of the stress-energy tensor - exactly the 
same reason that this constant isn't 'gauge' for the ordinary perfect fluid solutions with 
equation of state P = kp.) 

A metric is considered to be asymptotically flat [^,0 if in the vicinity of a spacelike 
hypersurface its components behave as 

aafe(x7r,t) 1 

gab ^ Vab + + (16) 

as r oo. (e > 0, rjab is the Minkowski metric, Oab is an arbitrary symmetric tensor, and 
X IS cl flat coordinate system at spacelike inflnity.) 

According to this deflnition, our metrics (|T5|) are asymptotically flat for k > 1/2 and are 
cosmological for k < 1/2. (In particular, for k = 1, f{v) = M and 2g{v) = Q^, the metric 
is just Reissner- Nordstrom). 

6 



When the imposed equation of state on the eigenvalues (|) and (J^) of the stress-energy 
tensor (RI) is P = kp"", the equation for the mass function is more comphcated : 



m" . , m' a 



k ij—S- (17) 



This has solution 

m{r,v) = fiv) + j dr[ g{v) - k {^Anf-^-^i^-'^) (Ig) 

Since a < 1, this mass function gives only cosmological metrics - the pressures are too small 
to make them asymptotically fiat. However if the dominant energy condition is not imposed 
but only the weak (or strong) one is, then a > 1 is possible, which will give asymptotically 
fiat metrics. 

We now consider in turn the asymptotically flat {k > 1/2) and the cosmological {k < 1/2) 
metrics. 

k > 1/2 : Since the dominant energy condition is satisfied only for g{v) > and g{v) < 0, 
if g{v) is zero initially, it must remain zero. On the other hand, if it is non-zero initially, 
it can decrease to zero. Therefore the asymptotically fiat metrics can be fiat in the v — >■ 
— oo limit only if g{v) = 0. But this is just the Vaidya case. However the following two 
parameter (A > 0, < -B < 1) family of metrics, (where we have set k = 1 for concreteness) , 
demonstrates an interesting feature: 

ds^ = -(1 - ^(^+^""^"^ + ^^-^f"^"^) dv^ + 2 dvdr + r-dn- (19) 

In the V — oo limit this metric has a naked singularity at r = 0. However, in the v +cx3 
limit it may have horizons depending on the relative values of A^B. Specifically, these 



horizons are given by r = A ± A'^ + B — 1. Therefore a black hole first forms at non-zero 
mass when A^ = 1 — B. This mass gap seperates a black hole from a naked singularity, and 
may be contrasted with the critical behaviour solutions where there is no mass gap 
between fiat space and a black hole metric. Similar results hold for other values of A; > 1/2. 
What is happening physically is that initially only the 'charge' term is present in the metric, 
and subsequently, infalling fluid reduces the 'charge' and adds a 'mass' term. 



These solutions also include metrics that give the evolution of a naked singularity at 
V = — oo into flat space at f = oo. This case occurs if we set A = 0,B = 1 in (|19|). 

Another feature of these solutions is that they give black holes with null fluid hair for 
1/2 < k < 1. An example of such a static solution results from the v —>■ oo limit of the 
metric 

ds^ = -(1 - ^l(i±^ + (^--^J^^^H ^ ^ 2 dvdr + r^dn\ (20) 
with B as the 'hair'. (We have put k = 3/4). 

k < 1/2 : The dominant energy condition is now satisfied for g{v) > and g{v) > 0. There- 
fore these metrics can be flat as f ^ — oo, and have black holes as f ^ oo. A specific two 
parameter {A,B > 0) example for k = 1/3 (electromagnetic radiation) is 

ds^ = -( 1 - C + A{l + t^nhv) _ Bjl + tanh.) ^ ^ dvdr + rHn\ (21) 

The apparent horizons in the f — > oo limit (with C = 0) are now given by the cubic 
equation (r — 2AY — (2i?)^r = 0. This equation always has a solution for the ranges of 
A and B allowed by the energy conditions. Thus, this solution describes the collapse of 
electromagnetic radiation (P = p/3) from flat space at f = — oo to a black hole aX v = oo. 



Another possibility occurs when C is negative in (0). Then, as for the > 1/2 case 
above, we have a naked singularity at f = — oo, which becomes a black hole at f = oo as 
the collapse proceeds. 

In the general case, if lim^^oo f{v) = A and lim^j^oo 5'('^) = B the radii of the apparent 
horizons in the f oo limit of the metric (pIS]) are given by 



2B 

{2k -^ryv^ 



2^ - .o. .^...-1 ' (22) 



which in general may have multiple solutions. 



In conclusion, we have given a new class of null fluid collapse solutions (|T5|) and (|T^) of 
the Einstein equations for the stress-energy tensor (|^). These include new asymptotically 
flat black hole solutions (a = 1,1/2 < k < 1) with multiple apparent horizons and hair. 
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The general metric depends on one parameter (k), and two arbitrary functions of v (modulo 
energy conditions). The long time limits of the asymptotically fiat solutions 'fall between' 
the Schwarzschild and Reissner-Nordstrom metrics in the sense that 1 < 2A; < 2 in (p!5|). 

Physically, the A; > 1/2 solutions describe the evolution of a naked singularity into the 
same, or into a black hole. The parameters in the solution (|T^ determine which of these 
possibilities occurs, and the black hole always forms at a finite non-zero mass. The A; < 1/2 
solutions describe the evolution of either fiat space or a naked singularity into a black hole 
in a cosmology. All of the new solutions support the cosmic censorship conjecture. 

Because of these physical properties, the exact solutions we have given are not of relevance 
for the collapse solutions that exibit critical behaviour The physical reason for this is 

that, although our stress-energy tensor (^) has non-zero components along both ingoing and 
outgoing null directions, there is energy flow only along one direction, (because, as noted 
above TabV°'v^ = 0). 

It should be possible to find more exact solutions by imposing other equations of state 
on the eigenvalues of the stress-energy tensor (^. 

This work was supported by the Natural Science and Engineering Research Council of 
Canada, and by NSF grant PHY 93-96246 to the Pennsylvania State University. 
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